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Abstract. We prove Holder continuity of the integrated density of states for 
the Fibonacci Hamiltonian for any positive coupling, and obtain the asymp- 
totics of the Holder exponents for large and small couplings. 



1. Introduction 

In this paper we are interested in regularity properties of the integrated density of 
states associated with the Fibonacci Hamiltonian. We will show that it is uniformly 
Holder continuous and provide explicit estimates for the Holder exponent in terms 
of the coupling constant. 

To motivate our study, let us consider an invertible ergodic transformation T of 
a probability measure space (f2, dfi) and a bounded measurable function / : fl — > R. 
One associates a family of discrete Schrodinger operators on the line as follows: For 
ui € 0, the potential : Z — > K is given by V w (n) = f(T n uj) and the operator 
in £ 2 (Z) acts as 

[H u <f>](n) = if>(n + 1) + 4>{n - 1) + V u (n)4>(n). 

An important quantity associated with such a family of operators, {i? w } w en, is 
given by the integrated density of states, which is defined as follows; compare 
[21 03] . Define the measure dN by 



(1) / g(X)dN(X) = J (So,g(H u )5o) dfi{v). 
The integrated density of states (IDS), N, is then given by 

(2) N(E) = J X(_ co , £] (A)diV(A). 

The terminology is explained by 

ATf J? , ,. #{eigenvalues of H Ui[1 , n] < E} 

(3) N(E) = hm — for /x-a.e. u € il, 

n— >oo Ti 

where H u ti >n ] denotes the restriction of H u to the interval [l,n\ with Dirichlet 
boundary conditions. It is a basic result that the IDS is always continuous see 
[2 6) for a very short proof that also works in higher dimensions. Craig and Simon, 
[TT] (see [12] for the multi-dimensional case), have shown that more is true. Using 
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the Thouless formula, they proved that N is log-Holder continuous, that is, there 
is a constant C such that 

(4) \N(E 1 ) - N(E 2 )\ < C (log \E, - E 2 \- 1 y 1 

for all Ei,E 2 with \E X -E 2 \<l/2, 

In this general setting, the bound Q is essentially optimal; see Craig jTUj and 
Gan and Kriiger |27j . However, for concrete models, one may hope to improve upon 
Q. Roughly speaking, one expects stronger regularity properties of the IDS the 
more random the stochastic process VL,(n) is. In the i.i.d. situation, it is always 
Holder continuous, that is, 

(5) \N(E 1 )-N(E 2 )\<C\E 1 -E 2 \i 

for some C < oo and 7 > 0, as shown by Le Page [37]. If the single-site distribu- 
tion is nice enough, Simon- Taylor |44j and Campanino-Klein [7] proved that N is 
C°°. We refer the reader to the survey article [13] by Simon, which describes the 
regularity results for the IDS that had been obtained by the mid-1980's. 

More recently, there has been renewed interest in the problem of proving regu- 
larity better than Q for some classes of ergodic Schrodinger operators with little 
randomness. This was initiated by Goldstein and Schlag, who proved Holder conti- 
nuity of the IDS, ([5]), for analytic quasi-periodic potentials in the regime of positive 
Lyapunov exponents [25]. This paper was followed by [H 03 [S] E2|] who improved the 
estimate in some cases or proved regularity results for different models (e.g., with 
T given by a skew-shift on T 2 ). Also see Hadj-Amor [SU] and Avila-Jitomirskaya 
[I] for results for analytic quasi-periodic potentials in the regime of zero Lyapunov 
exponents, and Schlag [32] for results for analytic quasi-periodic models at large 
coupling in two dimensions. 

Due to the typical presence of a dense set of gaps in the spectrum, one does 
not hope for more than Holder regularity for quasi-periodic (or, more generally, 
almost-periodic) models. 

Our objective here is to study the regularity properties of the IDS for a prominent 
quasi-periodic model that is not covered by the Goldstein-Schlag paper and its 
successors; the Fibonacci Hamiltonian, introduced independently by Kohmoto et 
al. [34] and Ostlund et al. [38]. It is given by = T = M/Z, Tx = x + a mod 1, 
where _ 

Vb - 1 

is the inverse of the golden mean, /1 is the Lebesgue measure on T, and f{uj) — 
^X[i-a,i) ( w ) f° r some A > 0. Thus, the potentials have the form 

(6) V u (n) = \X[i-a,i)( na + w mod !)• 

The associated operators {H u } ue Q form the family of Fibonacci Hamiltonians. 
This is the standard model of a one-dimensional quasicrystal. The spectrum of H u 
is easily seen to be independent of uj, and it will henceforth be denoted by £>. 
This set is known to be a Cantor set of zero Lebesgue measure [H] (and in fact of 
Hausdorff dimension strictly between and 1 [5]). Moreover, it is known that H u 
has purely singular continuous spectrum for every A and u>; see Damanik and Lenz 
[23] . Kotani [35], and Siito [45l [46]. The survey articles O QU [TBI El] contain 
information on the results obtained for this model and its generalizations. Here 
we are interested in the integrated density of states associates with the Fibonacci 
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Hamiltonian, which we will henceforth denote by N\ since its dependence on the 
coupling constant A will be of explicit interest. We mention in passing that dN\ is 
the equilibrium measure on £a in the sense of logarithmic potential theory. 

We first note a result that is essentially well-known, but which is stated for the 
sake of completeness. 

Theorem 1. For every A > 0, there are C\ < oo and 7a > such that 

{NxiEi) - N X (E 2 )\ < CxlEt - Etf* 
for every E\ , E 2 with \E\ — E 2 \ < 1 . 

Proof. It follows from the definition Q-Q that the integrated density of states 
is the distribution function of the /i-average of the spectral measures with respect 
to and 8q. It was shown in :52] that for each A, these spectral measures are 
uniformly Holder continuous with constants C\ < oo and 7a > that arc uniform in 
ui. That is, the /z-average of these measures will also be uniformly Holder continuous 
with the same pair of constants. □ 

One can infer explicit expressions for C\ and 7a from |22j . However, they are 
clearly far from optimal and hence we opted not to make them explicit. 

Our main goal is to identify the asymptotic behavior of the Holder exponent in 
the regimes of large and small coupling. As mentioned above, the tools used to 
establish Theorem [T] do not produce optimal results and hence are inadequate to 
identify the asymptotic behavior precisely. Therefore, different methods are needed 
in these asymptotic regimes, and we will indeed use different ones in either of these 
two cases. 

In the large coupling regime, we have the following: 
Theorem 2. (a) Suppose A > 4. Then, for every 

7 21og(2A + 22)' 

there is some 5 > such that the IDS associated with the family of Fibonacci 
Hamiltonians satisfies 

\NxiEt) - N X (E 2 )\ < \Ei-Ehp 
for every E\,E 2 with \E\ — E 2 \ < S. 
(b) Suppose A > 8. Then, for every 



7 > 



21og(i((A-4) + v /(A-4)2-l 2 )) 

and every < 5 < 1, there are Ei,E 2 with < \E\ — E 2 \ < S such that 

\NxiEt) - N X (E 2 )\ > \Ei-Erf. 

This shows in particular that the optimal Holder exponent is asymptotically 

31 2kjgA " m ^ ne l ar S e coupling regime. 

Theorem [2] is proved in Subsection |2.2| The proof is based on the self-similarity 
of the spectrum of H u . In particular, we do not use the Thouless formula and a 
Holder continuity result for the Lyapunov exponent, as was the case in many of 
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the works mentioned above. Thus, in a sense, we use a geometric, rather than 
an analytic, approach. Before turning to the proof, we first recall the canonical 
periodic approximants, which are obtained by replacing a by its continued fraction 
approximants. This enables us to describe the self-similarity of the spectrum that 
is crucial to our proof and it will also establish an explicit way to express N\(E) 
in terms of periodic spectra. 

In the small coupling regime, we have the following: 

Theorem 3. The integrated density of states N\(-) is Holder continuous with 
Holder exponent j\, where 7a — > ~ as A — > 0, and 7a < 5 for small A > 0. 
More precisely, 

(a) For any 7 G (0, there exists \q > such that for any A G (0, Ao), there 
exists S > such that 

\NxiEx) - N X (E 2 )\ < \Ex -£ 2 p 

for every Ex,E% with \E\ — E%\ < S. 

(b) For any sufficiently small A > 0, there exists 7 = 7(A) < | such that for 
every S > 0, there are E\, Ei with < \E\ — Ei\ < S and 

\Nx(Eh.) - N X (E 2 )\ > \Et-E2p. 

Theorem [3] is proved in Subsection |3.6| The proof uses the trace map formalism 
and the dynamical properties of the Fibonacci trace map studied previously in 
[31 [HI [HI ESI EE [23 EI]- In particular, a relation between the integrated density of 
states of the Fibonacci Hamiltonian and the measure of maximal entropy of the 
trace map was established in [21]. The proof combines this relation with Holder 
structures that appear due to hyperbolicity of the trace map in order to get explicit 
estimates on the Holder exponent. To show that the obtained asymptotics of the 
Holder exponent are optimal, we study the behavior of unstable multipliers of 
specific periodic points of the trace map. 

We conclude this introduction with some general remarks. In [21] we studied the 
scaling exponents associated with the measures dNx- We showed that there exists 
< A < 00 such that for A G (0, Ao), there is d\ G (0, 1) so that the density of 
states measure dN\ is of exact dimension d\, that is, for diV^-almost every £el, 
we have 

.. log 7V A (£ - £ ,E + e) 

hm = d\ . 

eio log e 

Moreover, 

limdx = 1. 

H0 

While at first sight this result and the question addressed in Theorem [3] are quite 
similar, note that in the latter result, one has to establish a uniform estimate in 
the vicinity of an arbitrarily chosen energy in whereas in the former result, 
one may exclude a set of zero dNx measure from the consideration. This crucial 
difference leads to different answers, as the different asymptotics show. Of course 
the scaling exponent of dNx at any energy E G Sa bounds the optimal global 
Holder exponent of Nx from above, and hence one gets a one-sided estimate in this 
way. More specifically, the scaling exponents are worse at gap boundaries of Ea 
and essentially determine the global Holder exponent of Nx , but these points form 
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a countable set and hence a set of measure zero with respect to dN\ since this 
measure is always continuous. 

It is interesting to compare the large coupling asymptotics of several A-dependent 
quantities. As we saw in Theorem [2j the optimal Holder exponent behaves asymp- 
totically like 1.5 • lo fo° A ■ On the other hand, it was shown in [17] that the 

Hausdorff dimension of £>, behaves asymptotically like 1.831 . . . • A ^ , and it 
was shown in [25] that a certain transport exponent, which measures the rate of 
wavepacket spreading in the time-dependent Schrodinger equation associated with 

the Fibonacci Hamitonian, behaves asymptotically like 2 • lo ^° A ^ . Similarly, in the 
weak coupling regime, we have that the Hausdorff dimension of T,\ strictly exceeds 
the dimension d\ of dN\ (a fact that was also proven in |21l ). which in turn strictly 
exceeds the optimal Holder exponent of N\ (since they have different asymptotic 
values). This shows that the strongly coupled and the weakly coupled Fibonacci 
Hamiltonian serve as a good source of examples demonstrating that certain quan- 
tities associated with a discrete Schrodinger operator need not be identical. More- 
over, the three different prefactors in the large coupling asymptotics (1.5, 1.831 . . ., 
and 2) correspond directly, and in a quite beautiful way, to the scaling properties 



exhibited by Ea (cf. Subsection 2.1 1 



Finally, determining the correct Holder exponent is one important ingredient in a 
recent study of the spacings of the zeros of a certain class of orthogonal polynomials 
(or, equivalently, the eigenvalues of H u \i ^ in the notation introduced above) by 
Kriiger and Simon [38j and hence our results feed into their theory. 

Acknowledgments. It is a pleasure to thank Mark Embree for help with symbolic 
calculations and Helge Kriiger and Barry Simon for useful comments. In particular 
we are grateful to Helge Kriiger for his question about the small coupling behavior 
of the Holder exponent, which prompted us to prove Theorem [3j 



2. The Large Coupling Regime 

2.1. Canonical Periodic Approximants and Scaling Properties. In this sub- 
section, we recall some known results for the Fibonacci Hamiltonian, its periodic 
approximants, and their spectra. The main tools we shall need in the sequel are 
summarized in Propositions |2.1f|2.3| below. 

Define the sequence {F k ) k >-i of Fibonacci numbers by 

F_ x = 0, F = 1, F k = F fe _! + F fe _2 for k > 1. 

For k > 1, define 

x k (E,X) = *tr M k (E,X), 

where 

-( s - v ->™ 1)*-*( E ~ v r (1) 1)- 

that is, the transfer matrix for phase uj = and energy E from the origin to the 
site F k . 

The matrices M k obey the recursion 



(7) 



M k (E, A) = M k _ 2 (E, X)M k _ 1 (E, A), 
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k + 2 

fc + 1 
k 

k-1 



A 2 



Figure 1. left: a type A band in a k ; right: a type B band in a k . 
and as a consequence, 

(8) x k+1 (E, A) = 2x k (E, X)x k -i{E, A) - x k - 2 {E, A). 

With the definitions above, ^ and Q hold for k > 3. If we define 

M^(E,X) = ( J ~M and M (£,A) = ( f ~J 

these recursions extend to all fc > I. 

The so-called trace map relation, pj), yields the invariant 

(9) 

z fe+ i(£, A) 2 + z fe (£, A) 2 + x k -i(E, A) 2 - 2a: fe+1 (£;, X)xh(E, X)x k -i(E, X) = I + 

The identities Q-© were proved by Siito in |45j . 

For fixed A, define (leaving the dependence on A implicit) 

a k = {EeR:\x k (E,X)\<l}. 

The set a k is actually equal to the spectrum of the Schrodinger operator H whose 
potential V k results from Kj=o in ^ by replacing a by Fk~i/F k (see [IS]). Hence, 
Vk is i^fe-periodic, a k C R, and it consists of F k bands (closed intervals). 

Next, we recall some results of Damanik and Tcheremchantsev 24;, Killip, Kise- 
lev, and Last [33], and Raymond 40J. From now on, we shall always assume A > 4 
since we will make critical use of the fact that in this case, it follows from the 
invariant, Q, that three consecutive half-traces cannot simultaneously be bounded 
in absolute value by 1: 

(10) VA>4, VE,k : m&x{\x k {E,X)\,\x k+1 (E,X)\,\x k+2 {E,X)\} > 1. 

Following [33], we call a band I k C a k a type A band if I k C <y k ~i (and hence 
Ik H (cfc+i U <J k --2) = 0). A band I k C cr k is called a type B band if I k C cr k -% (and 
therefore I k n cx k -x — 0); compare Figure 1. 

By definition of M_i and M , cr_i = M, cr = [-2, 2], and in = [A - 2, A + 2]. 
Hence, <7o consists of a single band of type A, and o\ consists of a single band of 
type B. 



From (10 1, one gets the following result on the relative position of bands on 
successive levels (cf. [331 Lemma 5.3] and [101 Lemma 6.1]): 

Proposition 2.1. For every A > 4 and every k > 0, 

(a) Every type A band I k C o~ k contains exactly one type B band 1^+2 C o- k+2 > 
and no other bands from a k+ i, o~ k +2- 
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(b) Every type B band Ik C ut contains exactly one type A band Ik+i C crfc+i 
and two type B bands from o~k+2, positioned around Ik+i- 

The Lebesgue measure of each individual band admits the following geometric 
(in k) lower bound: 

Proposition 2.2. For A > 4, k > 3, and every band Ik of Ok, we have 

4 



where 



Ik > 



denotes Lebesgue measure. 



(2A + 22) 2fc / 3 ' 



Proof. This follows from [2H Lemma 3.5], noting (10) and the fact that on each 
band Ik of o- k , we have dE = 2^ 



□ 



In the large coupling regime, this result is asymptotically optimal in the following 
sense: 

Proposition 2.3. For A > 8, k > 4 and k = 1 mod 3, there exists a type B band 
Ik of o~k with 

4 

4 < 



(i ((A-4) + V(A-4)*-12)) 



2fc/3 ' 



Proof. This follows from [331 Lemma 5.5], Proposition |2.1| above (start with the 
type B band of o\ and then cycle periodically through B, A, and "not contained"), 
and again the fact that on each band Ik of <7fc, we have dE = 2. □ 

Let us study the scaling properties of the sets <Jh in more detail. In particular, 
we want to show that, when n > k, each band of o~k contains either F n -k or F n —k—2 
bands of a n . 

Recall that ctq consists of a single band of type A, and o\ consists of a single 
band of type B. A repeated application of Proposition |2 . 1 1 yields the entries in the 
following table: 



k 


# of bands in of type A 


# of bands in o~k of type B 





1 





1 





1 


2 


1 


1 


3 


1 


2 


4 


2 


3 


5 


3 


5 


6 


5 


8 









We observe the following: 



Here we caution the reader that the literature on the Fibonacci Hamiltonian is not consistent 
in the sense that some papers use half-traces, as we do here, and other papers use traces instead. 
The papers 24 33 40 , on which much of the present subsection is based, belong to the second 
group and hence their results need to be slightly reformulated when given here. 
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Lemma 2.4. For A > 4 and k > 2, we have 

# of type A bands in o~k = -Ffc-2, 

# of type B bands in o~k = i*fc-i- 

In particular, for k > 0, every band of o~k is either of type A or of type B. 
Proof. This is a straightforward induction. □ 
A similar straightforward application of PropositionpO]gives the following result. 



Lemma 2.5. Suppose A > 4, k > 0, and n > k. 

(a) Every type A band of contains F n —k-2 bands of o~ n . 

(b) Every type B band of contains f n — fc bands of o~ n . 

2.2. Proof of Theorem [2| In this subsection we prove Theorem [2] O ur mam 



tools will be Proposition 2.2 Lemmas 2.4 and 2.5 and a formula, given in (11) 



below, connecting the IDS of the Fibonacci Hamiltonian and the band structure of 
the periodic spectra introduced the previous subsection. 

Let us first note that in the Fibonacci case, as a consequence of results of Hof for 
uniquely ergodic models, [31] . the convergence in Q takes place for every, rather 
than almost every, uj e il. Moreover, Dirichlet boundary conditions can be replaced 
by other boundary conditions, such as Neumann or periodic boundary conditions. 
If we choose periodic boundary conditions, and consider u = and convergence 
only along the subsequence (-F n )n>i, we obtain the following formula (which had 
already been noted by Raymond [40 J, 

(11) N X (E)= hm #{bandsofq»<^ 

n— >oo F n 

since each band of a n contains exactly one eigenvalue of the operator H u= q, re- 
stricted to [l,F n ] with periodic boundary conditions. 

Proof of Theorem^ (a) Given E 1 < E 2 (with E 2 -E 1 small; less than 4/(2A+22) 2 , 
say), we want to estimate N\(E2) 

(12) N\(E 2 ) — N\(Ex) = lim 



say), we want to estimate N\(E 2 ) — N\(E{) from above. It follows from (111 that 

bands of a n contained in [Si,^]} 

Thus, we need to estimate from above the number of bands of a n that are contained 
in [Ei, Ea]. 

Let k > 3 be the integer with 

4 4 
< E 2 - Ei < 



(2A + 22) 2 ( fc + 1 )/ 3 " " 1 (2A + 22) 2fc / 3 ' 

Consider, for n > k, the bands of a n that are contained in [Ei,E 2 j. By the results 
of the previous section, each of these bands is contained in a band of or in a 
band of o~k-i- By the definition of k and Proposition 2.2[ at most two bands from 
(jfe U (Jfe_i can occur as associated bands. Thus, Lemmas 
most 2F n _(u_i\ bands of o~ n can be contained in [Ei,E 2 ] 
Define by 

(fc-3)log(a- 1 ) 



2.4 and 2.5 imply that at 



7fe = 2 



(fc + l)log(2A + 22)-log4 
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Then, we infer from (12 1 that 

„ ,. #{bands of er„ contained in [E 1 ,E 2 ]} 
N\{E 2 ) - Nx(Ei) = hm 

n— >oo r n 
< hm 2F - (fc - 1) 



< lim 



ra-(fc-3) 



n— ¥00 F n 
„ k-3 



Now, given 



,(2A + 22) 2 ( fc + 1 )/ 3 
< (E 2 -E 1 y k . 



o< 7 < 31og ^ 



21og(2A + 22)' 
choose ko > 3 such that 7fc > 7 for every k > ko. Let 

<5= (2A + 2V/ 3e(0 ' 1) - 
We obtain that for every E\,E 2 with \E\ — E 2 \ < 5, 

\NxiEj.) - N X (E 2 )\ < \E, - E 2 \i» < \Ex - E 2 \\ 

where k > k is the integer associated with \E\ — E 2 \. 
(b) Suppose A > 8 and let 



76 

and < 8 < 1 be given. 
Choose k such that 



Slog^" 1 ) 



21og(§ ((A-4) + v /(A-4)2-12 



2fc /3 



< <S 



(5 ((A-4) + V(A-4)2-12)) 
and such that 

(fc + l)loga _1 

7fc := 7 — 7 — < 7 

ffc log (I ((A - 4) + v/(A-4)2-12)) - log4 

for every k> k$. 

Next, choose fc > max{/co,4} with fc = 1 mod 3. By Proposition 2.3 
exists a type B band Ik of <Jk with 

4 

/* < 



(l((A-4) + v /(A-4)2-l 2 )) 2fc/3 



Denote the endpoints of Ik by Ei,E 2 , that is, I k = [E\, E 2 \. 
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By ( 12 ), Lemma 2.5 the definition of jk, the choice of Ik = [E±, E 2 ], the choice 
of k, and the fact that < E 2 — E\ < S < 1, we find that 

7^{bands of er„ contained in [Ei,E 2 ]} 



N X (E 2 ) - Nx(Et) = lim 



= lim 

n— >oo 

= a k+1 



F-n-k-l 



> (E 2 -E 1 y" 

> {e 2 -e 1 )\ 



12 



2fc/3 



This completes the proof. 



3. The Small Coupling Regime 

Here we will use the relation between the IDS for the Fibonacci Hamiltonian 
and the measure of maximal entropy for the so called Trace Map associated with 
the discrete Schrodinger operator with Fibonacci potential. 

3.1. IDS for the Free Laplacian. It is well known that So = [—2,2] and 

'0 E<-2 
arccos(-f) -2 < E < 2 
1 E > 2. 



(13) 



N (E) 



In particular, Nq is Holder continuous with Holder exponent 1/2. We will need a 
somewhat more detailed description of the continuity properties of A^o (which also 
follows directly from the explicit form (13)). 

Lemma 3.1. For any e > 0, there exists C > such that the following hold. 

(a) IfE 1 ,E 2 e [-2 + e,2 - e], then 

\N (E 2 ) - N (E 1 )\ <C\E 2 -E 1 \. 

(b) IfE u E 2 e [-2,2], then 

\N Q (E 2 ) - AT (£i)| < C\E 2 - E1I 1 / 2 . 

(c) IfE 1 ,E 2 e [-2, -2 + e) and -2 < E 1 < E 2 , then 

C 



\N (E 2 ) - NoiEt)] < 



|2 + ia|i/2 

Similarly, if Ex, E 2 G (2 — e, 2] and Ei < E 2 < 2, then 

C 



E 2 -E l \ 



\N Q {E 2 ) - N {E X )\ < 



\E2-EA 



|2-E 2 |V2' 

(d) There exists Cq > such that for any E 6 (—2, 2), one has 
\N {-2)-N {E)\>C \E + 2\ 1 / 2 and \N {2) - N Q {E)\ > C \E - 2] 1 ' 2 
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3.2. The Trace Map. The recursion ^ gives rise to a fundamental connection 
between the spectral properties of the Fibonacci Hamiltonian and the dynamics of 
the trace map 

T : E 3 ->• E 3 , T(x, y, z) = (2xy - z, x, y). 

The function G(x, y, z) = x 2 + y 2 + z 2 — 2xyz — 1 is invarianirl under the action of 
T (cf. ([9])), and hence T preserves the family of cubic surface^] 

f A 2 

S x = < (x, y, z) e E 3 : x 2 + y 2 + z 2 - 2xyz = 1 + — 

It is therefore natural to consider the restriction T\ of the trace map T to the 
invariant surface S\. That is, Tx : S\ — > Sx, Tx = T\g x . We denote by Ax the 
set of points in Sx whose full orbits under Tx are bounded (it is known that Aa is 
equal to the non- wandering set of T\). 

3.3. Hyperbolicity of the Trace Map. Recall that an invariant closed set A of 
a diffeomorphism / : M — ¥ M is hyperbolic if there exists a splitting of the tangent 
space T X M = E^®E^ at every point a: € A such that this splitting is invariant under 
Df, the differential Df exponentially contracts vectors from the stable subspaces 
{E%.}, and the differential of the inverse, Df -1 , exponentially contracts vectors from 
the unstable subspaces {E%}. A hyperbolic set A of a diffeomorphism / : M — > M 
is locally maximal if there exists a neighborhood U of A such that 

A = f| /"([/). 

It is known that for A > 0, Aa is a locally maximal hyperbolic set of T\ : S\ — > S\; 

see Banana. 

3.4. Properties of the Trace Map for A = 0. The surface 

§ = S n {{x, y, z) e E 3 : |z| < 1, \y\ < 1, \z\ < 1} 

is homeomorphic to S" 2 , invariant under T, smooth everywhere except at the four 
points P 1 = (1,1,1), P 2 = (-1,-1,1), ^3 = (1,-1,-1), and P 4 = (-1,1,-1), 
where S has conic singularities, and the trace map T restricted to S is a factor of 
the hyperbolic automorphism of T 2 = E 2 /Z 2 given by 

A(d,<p) = {d + ip,6) (mod 1). 

The semi-conjugacy is given by the map 

(14) F : (6,tp) 1 y (cos2tt(6I + (p), cos 2tt6>, cos 2ntp). 

The map A is hyperbolic, and is given by the matrix A — ^ j , which has 
eigenvalues 

l + \/5 , _j 1-V5 

a = and — a = . 

P 2 '2 

A Markov partition for the map A : T 2 —> T 2 is shown in Figure [2] Its image 

under the map F : T 2 —> S is a Markov partition for the pseudo-Anosov map 

T : § -> S. 

2 It is usually called the Fricke-Vogt invariant. 
3 The surface So is known as Cay ley cubic. 
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(0,1) 



(1,1) 



(0,0) 




-0.5 



Figure 2. The semi-conjugacy F between the linear map A and 
the trace map T on the central part S of the Cayley cubic. 

3.5. Spectrum and Trace Map. Denote by l\ the line 



It is easy to check that l\ C S\. An energy E £ K belongs to the spectrum 
T,\ of the Fibonacci Hamiltonian if and only if the positive semiorbit of the point 
(^T^'lT'l) un der iterates of the trace map T is bounded; see [45]. Moreover, 
the stable manifolds of points in A\ intersect the line £\ transversally if A > is 
sufficiently small [IB] or if A > 16 [5]. 

Let us denote L\ : K — > t\,L\{E) = (^p^f^l)- It is affine and contracts 
distances by the multiplicative factor ^ . 

Define the map ^a : S A -> [0, 2] by 



It turns out that there is a direct relation between the map \&a and dynamical 
structures of the trace map T\. The following statement is implicitly contained in 
gH Claim 3.2]. 

Proposition 3.2. There exists Xq > such that the following holds. Take any 
As [0,Ao) andx\,X2 £ ^a- Consider the stable manifolds W s {x{) andW s (x2), and 
take some points pi = W s (x\) n i\ and P2 = W s (x2) fl l\. When A changes, there 
are unique continuations of the points Xi,X2 £ denote them by Xi(A), 0:2 (A). 
The continuations of the intersections Pi(X) = W s (xi(X)) D l\, i = 1,2, are also 
well defined, and the value of the difference N\(L' X (p2(A)) — N\(L~^ (j>i(A)) is 
independent of X £ [0, Ao). 

Notice that Proposition |3.2| gives a dynamical description of the map ^a- In 
[2Tj this description was used to establish a relation between the IDS N\ and the 
measure of maximal entropy for T\\ x . 




(15) 



=V N x {x) = N {y). 
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3.6. Holder Continuity of the IDS in the Small Coupling Regime. Here 
we prove Theorem [3j 

Proof of Theorem^ Choose sufficiently small neighborhoods U(Pi) of the singu- 
larities {Pi, P2, P3, P4} of the Cayley cubic. Let U*{Pi) C U(Pi) be an essen- 
tially smaller neighborhood of the singularity P t . Set U — U J= i 2 3 4 ^(-^) anc ^ 
U* = U- =1 2 3 4 U*(P l ). Take any E 1: E 2 € Ex, and denote by b the interval on l\ 
between the points La (Pi) and L\{E-i). 
Notice that 

I^a(Pi) - N X (E 2 )\ = |JV (*a(-Bi)) - N (V X (E 2 ))\. 

Let a be the interval on Iq between the points Lo(5'a(Pi)) and Lq(^\{E<x)). Wc 
will consider separately three cases: 

(i) a and b are away from the neighborhoods J7*(Pj) of the singularities. 

(ii) a and b are in a neighborhood U{Pj) of a singularity (that could be either 
Pi = (1, 1, 1) or P2 = (—1, —1, 1)), and one of the edges of a is on a local 
strong stable manifold of the singularity (which implies that one of the 
edges of b is on the local strong stable manifold of a periodic orbit of period 
2 or 6) . This is equivalent to the case when Pi or P2 is equal to min T,\ or 
max Ea- 

(iii) a and b are in a neighborhood P(Pj) of a singularity (that could be either 
Pi — (1, 1, 1) or P 2 = (—1, —1, 1)), and none of the edges of a is on a local 
strong stable manifold of the singularity. 

Note that if |Pj — P 2 | is sufficiently small (depending on the choice of U and 
U*), then exactly one of the cases (i)-(iii) applies. 

Consider the case (i). Let us iterate a and b until they grow up to the length of 
order one. To simplify the estimates, let us introduce the following notation: 

x~ c y C- X |y| < \x\ < C\y\. 

There exists a constant C > that is independent of a, b such that for some 
M 6 N, we have \T M (a)\ ~ c |P M (&)| ~c L Let us split the sequence 
of iterates {a,T(a),T 2 (a), . . . ,T M (a)} into finite intervals {a, T(a), . . . , T kl (a)}, 
{P fel+1 (a), . . .,T fc2 (a)}, . . . (and similarly for {b, T(b), T 2 (6), . . . ,T M (b)}) in such a 
way that for each ki, one of the following cases holds: 

(a) {T fc - 1+1 (a), . . . ,T fc »(a)} as well as {P fc - 1+1 (&), . . .,T ki (b)} are away from 

U*, and moreover, ^Ciff), > 2 and > 2 ' 

(b) {T fc *- 1+1 (a), . . . ,P fc *(a)} as well as {T fc - 1+1 (6), . . . ,T fe *(&)} are inside of U. 

In the case (a), since for small A > 0, the maps Pols\c/* and P\|s A \(7* are en- 
close, for some a A < 1 with a x -» 1 as A -» 0, we have J^ziT^j > ( \T k ^(b)\ ) ■ 
In the case (b), by [20, Proposition 3.15], we have ' k _ " > Cfj, 2 , where (1q 
is an unstable multiplier at a singularity, and at the same time we have ^j, ki _\ ^ < 

C/i^ 1 where /i\ is an unstable multiplier of a periodic orbit Pa, and which 

obeys ^a — > A*o as A — > 0. Choose y\ < |, 7x — > \ as A — > 0, in such a way that 
/i^ 7a < /iQ- Then, (since in this case, the smallness of U* relative to U guarantees 



that ki - ki-i is large enough) we have j^^Mr, > ( J^t-i 



ML 
(6)1 



7A 
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Now we have 

1 ~ c \T M {a) 



\T*{a)\ \T"'(a)\ \T k ™(a)\ 



and 



Therefore, 



l~c \T M (b)\ = \b\ 



\a\ \T*(a)\ \T^(a)\ 

|T fcl (b)| \T k ?(b)\ \T k ™(b)\ 
\b\ ■ |T*x( 6 )| ••••• |r* — x( 6 )|- 



as) h= |rM( : )l < 

V 7 1 1 jjm |T^(g)| " 



Lj=l | T fcj-l( a ) 

< ~, ~u = ~r ^TtT < C 1+ ^\b\^ = C'\bp\ 

~ / r rm \T k Hb)\ \ 7A / |T^(b)| \ 7A - 11 11 

11.7=1 |T*J-i(6)| J { [51 J 



Since |6| = 4||-Ei - ^| and |a| = -^|* A (-Ei) - %\(Bh)\, using Lemma 3.1 (a), we 
find 

\N x (Ei) - 7Va(^ 2 )| = \N (V x (Ei)) - N (* X {E 2 ))\ < 

< C|* A (Bi) - *a(-B 2 )| < CV2|a| < CCv^lftl 7 * < C'l^i - E, 



Now let us consider the case (ii). In this case, one of the edges of a and b will 
never leave the neighborhood U* . Let us iterate a and b sufficiently many times to 
have 

(17) \T M (a)\ ~ c \T M (b)\ ~ c 1. 

Consider the coordinate system in a neighborhood of the singularity that rectifies 
all the invariant manifolds (see [HI Section 4] or (20j Subsection 3.2]). The lines 
£q and £\ are transversal to the central-stable manifolds, hence we can apply [2TJ1 
Proposition 3.15] with some bounded (independent of a and b) number k* from [201 
Lemma 3.16]. Since fx^ x < fJ,Q, we have 

|o| < Cmo M < C>a 27aM < C'|6| 27A . 
Using Lemma [3.1 |(b), we find 



\N x (Eh) - N X (E 2 )\ = \N (p\(.Ei)) - N (*x(E 2 ))\ < 

< C\*x(Ei) - *x(E 2 )\ 1/2 < C'\Ei - E 2 p\ 

Finally, let us consider the case (hi). Suppose that E\,E 2 G £a are such that 
Ei < E 2 < maxSx and ^x(Ei) is close to 2 (the case when ^x(Ei) is close to —2 
is similar). Then, a C £q> b C lx are in the neighborhood U(P\) of P\ = (1, 1, 1), 
and the distance from a to the central-stable manifold is of order fiQ S , where s is 
a number of iterates needed for a to leave U(P±); see [20 , Proposition 3.14]. This 
implies that 

|2-*a(£ 2 )|~ c 

Consider the arcs a — T s (a) and b = T s (b). These arcs are away from U, and 
therefore from (16) we deduce that \a\ < C|&| 7A . On the other hand, we have 

\a\ ~ c (i s \ a \, \b\ ~ c l*xK 
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and therefore 

\a\ < C^ s \a\ < CVo S H 7A < C 2 +T>o>riM 7A = C 2+7 Vo * (^%) H 7A < 



Now, using Lemma 3.1 (c), we find 



\N x (E t ) - N X (E 2 )\ = |7V (* A (^)) - JV (* A (^))I < 

^ 10 t^m^ '^^) - *^ i?2 )l ^ ^S-^o § H 7A < <?"|£i - BaP- 

Finally, notice that in all cases (i), (ii), (iii), the constant in the inequality can 
be set to 1 if one takes a slightly smaller Holder exponent and sufficiently small 
\E\ — E 2 \. This finishes the proof of part (a) of Theorem [3] 

In order to show part (b) of Theorem [3j consider the periodic points of period 2 
that are born out of the singularity (1, 1, 1). The strong stable manifolds of these 
points correspond to the boundaries of the gaps in the spectrum. These periodic 
points form a curve 

Per 2 = \ (x, y, z) : x g ( — oo, - ) U ( -, oo ) , y = X z = x, 



2 J \2 J 2x-l 

'<6 3 -2> 

For the map T 2 , these points are fixed points, and DT 2 (1, 1, 1) = I 2 2 — 1 I , 



1 



with eigenvalues 1 and 7± 3\/5 



Lemma 3.3. For A > sufficiently small, the largest eigenvalue of DT 2 at the 
periodic point of period 2 near the singularity (1,1,1) is strictly larger than the 
largest eigenvalue of DT 2 (1,1,1). 

Proof. Take a periodic point (x, 2x _i > € Per 2 . We have 



DT 2 [ x, 



2x-l 



1 2x(4x-l) 4x-l 2x 

2x-l (2x~l) 2 1-2.T 

^ X 2x~=l _ 
1 



V 



and the largest eigenvalue is 

. nt r \ 1 - 2x + 8x 2 + V-3 + 12a; + Ax 2 - 32x 3 + 64x 4 
{X) = ' 

Now we have 

^(*)U =1 = 0j ^( X )\ x=1 = 

_ 8(-3(l + 3^) + 2(9 + 18^5 + 2(-27 + 3\/5 + 2(55 - 12%/5 + 4(-l + 6\/5))))) 

135V5 

= 16.247987... > 

□ 
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In other words, Lemma |3.3| claims that fi\ > fj,Q if A > is small enough. Fix 
a small A > and take any 7 6 (0, |) such that /i^ 7 > /in. We claim that if 
E2 = maxEj, E\ £ Y,\, and \E\ — L2I is sufficiently small, then 

\N x (Ex) - N x {Ei)\ > \Bh.-E^. 

Indeed, ^(£'2) = 2, and the interval a C £q between the points Lr^vE^L^)) and 
Lo{^f x (Ei)) has one of its end points at the singularity Pi — (1,1,1). Consider 
also the interval b C £ x between La (£2) an d L X (E\). As in the case (ii) above, 



consider M iterates of a and 6, where M is such that (171 holds. Then, due to 
Proposition 3.14], we have 

Using Lemma |3.1| (d), we find 

\N x (Ek) - N X (E 2 )\ = |^ (*a(Li)) - N Q (^ X (E 2 ))\ > 



> C |*a(Li) - *a(L 2 )| 1/2 - 2 1 / 4 C |a| 1 /2 > C ' 




Mo 



>C" ^ \b\i>C" ^ IL2-L1I 7 > IL2-L1I 7 , 
Mo y V ^° / 

provided M is large enough (i.e., if |a| and |6| are small enough, or, equivalently, if 
\E2 — Ei I is small enough). This completes the proof of Theorem [3] □ 

Notice that as a byproduct of this proof, we also get the following statement: 





Proposition 3.4. The map : Ea —> [0,2] given by (151 is Holder continuous 
with a Holder exponent 7a that obeys j\ — > - as A — > 0. 

Remark 3.5. In terms of the dynamics of the trace map, there exists Xq > such 
that for any A € (0, Ao), the semiconjugacy $a : ^x — > &x £a|a a = ^ols 
is Holder continuous with a Holder exponent 7a such that 7a — > | as X — > 0. Since 
we are not using this statement here, we do not elaborate on it. 

This is related to the following classical fact from hyperbolic dynamics. Suppose 
Af is a compact locally maximal hyperbolic set of a surface diffeomorphism f : 
M — > M , g is C 1 -close to f, and A g is a continuation of Af. Then there is a 
continuous conjugacy h : Af —¥ A g , ho f = goh, and h has to be Holder continuous 
with Holder exponent close to one (see |32j . |39j). In our case, the Holder exponent 
of the conjugacy $a is close to 1/2, not to one, due to essentially different behavior 
ofTo\s andTx\h x near the singularities of the Cayley cubic. 
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